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ABSTRACT 

v  The  estimation  of  arrival  direction  is  an  important  task  in  signal 
processing  and  has  recently  received  considerable  attention  in  the  literature. 
In  this  paper,  the  authors  proposed  a  method  to  estimate  the  direction  of 
arrival  and  proved  the  strong  consistency  of  the  estimates  for  both  cases  in 


f 


presence  of  white  noise  and  colored  noise. 

u  , 

i  f  t  s  r,  ’j-  *  >  <  -■ 


AMS  1980  Subject  classifications:  Primary  62F12;  Secondary  62H12. 


Keywords  and  phrases:  Direction  of  arrival,  signal  processing. 


1 .  INTRODUCTION 


Since  the  work  of  Schmidt  (1981)  and  that  of  Bienvenu  (1979),  which 
in  turn  were  extensions  of  Pisarenko  (1973),  the  eigenstructure  methods 
for  direction  of  arrival  (DOA)  have  been  developed  rapidly  in  the  past  few 
years,  and  have  attracted  considerable  interest.  When  the  additive  sensor 
noise  is  spatially  white,  Wax,  Shan  and  Kailath  (1984)  proposed  a  method 
for  estimating  the  DOA.  This  method  is  based  on  the  fact  that  the  DOA 
vectors  are  orthogonal  to  those  eigenvectors  of  the  true  covariance  matrix 
of  observations  associated  with  the  smallest  eigenvalue.  In  some  cases, 
the  noise  is  not  spatially  white  and  its  covariance  is  unknown  and  in  this 
case  the  algorithm  of  Wax,  Shan  and  Kailath  is  no  longer  appliable.  In  these 
cases,  Paulraj  and  Kailath  (1986)  proposed  a  method  to  estimate  the  DOA  based 
on  the  difference  of  two  covariance  matrices.  Their  method  relies  on  the  fact 
that  the  DOA  vectors  are  orthogonal  to  the  eigenvectors  of  the  difference 
matrix  associated  with  the  zero  eigenvalue.  Both  methods  of  Shan -Wax -Kailath 
and  Paulraj -Kailath  are  based  upon  finding  the  infimum  of  a  Hermitian  form 
with  constrained  variables. 

However,  though  simulation  results  strongly  supported  the  above  two  methods 
for  estimating  the  DOA,  it  is  not  an  easy  task  to  find  the  solutions  for  the 

infimum  of  the  constrained  Hermitian  form.  In  the  present  paper,  we  investigate 

the  estimation  of  DOA  for  both  cases  where  the  noise  is  white  or  colored.  In  the 

argorithm  for  estimating  the  DOA,  we  only  need  to  solve  a  polynomial  equation 

whose  degree  is  just  the  number  of  signals.  Also,  we  shall  prove  that  this 
estimate  is  strongly  consistent  under  minor  moment  restrictions.  In  another 
paper  (in  preparation)  we  shall  investigate  the  asymptotic  normality  of  these 


estimates. 

The  organization  of  this  paper  is  as  follows:  In  Section  2,  we  shall 
describe  the  algorithm  for  estimating  DOA  when  the  noise  is  spatially  white 
and  prove  the  strong  consistency  of  these  estimates.  In  Section  3,  we  shall 
briefly  describe  the  procedure  for  finding  the  estimate  of  number  of  signals 
by  using  information  theoretical  criteria  and  the  estimate  of  DOA  by  the 
proposed  method  when  the  noise  is  colored.  We  only  point  out  that  these 
estimates  are  also  strongly  consistent  and  omit  the  details,  because  the 
proofs  are  almost  the  same  as  the  proof  for  the  strong  consistency  of  signal 
number  estimate  (see  Zhao,  Krishnaiah  and  Bai  (1986  a,b)  and  the  proof  given 
in  Section  2  for  the  strong  consistency  of  estimtes  of  DOA. 
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2.  ESTIMATE  OF  DOA  IN  THE  PRESENCE  OF  SPATIALLY  WHITE  NOISE 


Consider  the  model 


x(t)  *  As(t)  +  n(t). 


1 , 2, . . . ,N, 


(2.1) 


where  x(t):  pxl,  the  observations  received  by  p  sensors,  s(t):q*l,  the  signal 


vector  emitted  by  q  sources,  q  <  p,  n(t)  is  the  white  noise  vector,  A  *  (a^ . a^) 


,  "J<Vk 

and  a^=(l,e  ,...,e 


-jaj0(p-l)Tk  T 

)  ,  called  the  direction-frequency  vector 


associated  with  the  k1"*1  signal  j  =  /-T,  the  center  frequency  of  signals  and 


=  —  sin  0^,  A  the  spacing  between  sensors,  c  the  speed  of  propagation  and 

th 


0^  the  direction  of  k  signal.  Since  u)q  is  known,  we  can  assume  Wq  in  the  sequel 
It  is  usual  to  assume  that 

(i)  (s(t)}  are  independent  and  identically  distributed  (i.i.d.),  {n(tT}^ 

are  i.i.d.,  and  independent  of  {s(t)>  I 

(ii)  Es(t)  *  0, En(t)  '=  0,  Es(t)s*(t)  =  T  >  0, 


(2.2) 


2  2 

En(t)n*(t)  *  o  I  with  a  unknown, 
P 


(III)  Tjc's  are  distinct. 


where  *  denotes  complex  conjugate  transpose. 


Under  the  model  (2.1),  our  problem  is  to  find  an  estimate  of  t^'s  based 


on  the  sample  covariance  matrix 


N 


E  *  r  I  x(i)x*(i). 


N 


i-1 


The  covariance  matrix  of  x(t)  is  given  by 


ATA*  +  o  I 


Denote  by  A^  >_  X^  ^  X  and  ^2  —  "'  —  *^p  t^ie  e^8envalues  °f  2  and  E 

respectively.  Also,  let  e^,...,e  and  u^,...,u  denote  the  eigenvectors 
associated  with-  these  A's  and  6's  respectively.  Without  loss  of  generality,  we 
assume  that  u's  are  of  unit  length  and  orthogonal  of  each  other,  and  the  same 
is  true  for  e's.  If  the  number  of  sources  is  q,  then  we  have 

2 

X,  >  A_  >  ...  >  A  >  A  =  ...  =  A  (=0  ). 

1  —  2  —  —  q  q+1  p 

The  key  steps  of  Wax-  -  Shan  -  Kailath  algorithm  are  as  follows.  First, 

determine  the  number  of  sources  q.  Next,  find  the  so-called  noise  subspace  as 

the  span  of  the  eigenvectors  corresponding  to  the  minimal  (noise)  eigenvalues 
2 

a  of  E.  The  subspace  spanned  by  the  direction  vectors  of  the  impinging  signal 
wavefronts,  which  is  called  signal  subspace,  can  be  obtained  as  the  orthogonal 
complement  of  the  noise  subspace.  For  determining  the  DOA's,  they  plotted  the 
inverse  Hermitian  form  that  measures  the  orthogonality  between  the  direction 
vectors  and  the  noise  subspace,  i.e., 


(2.3) 


where  afl  =  (i,e  ^®,...,e  and  is  an  p*(p-q)  matrix  whose  columns 


are  the  eigenvectors  associated  with  the  minimum  eigenvalues  of  £.  They  pointed 

(Z)  (E) 

out  that,  "Ideally,  aQlE  ,  for  9=9,,  and  hence  H„  should  become  very  large 

n  k  y 

at  these  9^,  enabling  us  to  pick  out  the  source  directions."  In  other  words, 

(E)-l 

they  might  extract  these  9,  's  by  seeking  for  the  extreme  points  of  Ha  ,  a 

1C  o 

.40 

polynomial  of  e  J  with  degree  2(p-l). 


But  there  are  two  problems:  (1)  We  do  not  know  the  number  of  the  extreme  points. 
(2)  No  method  is  proposed  to  extract  the  desired  q  6.  's  from  these  extreme  points. 
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Now  we  introduce  a  new  method  as  follows: 

Since  the  information  theoretic  criterion  (ITC)  gives  strongly  consistent 
estimate  of  the  number  of  signals,  we  can  assume  that  q  is  known  throughout 
this  section,  (refer  to  Zhao,  Krishnaiah  and  Bai  [1986a]). 

Write 


W 


N 


(2.4) 


By  a  knowledge  of  linear  algebra,  there  exists  an  unitary  matrix  0^: (p-q)x(p-q) , 
such  that 


W0.=(u,,,...,u)  =  (u.,),.  ,,, 

n  N  -q+1  -p  lk  l£i£p,  q+lfk<p 


with  =  0,  for  k  =  q+1 . p-1  and  i  =  k+l,...,p, 


(2.5) 


uii  >  0,  for  i  =  q+1 , , 


Also,  if  all  u. .  >0,  then  0..  is  uniquely  determined, 
li  N 


Let  =  p^  exp(j?jt),  k  =  l,2,...,q  be  roots  of 


q+1 


B(z)  -  ^ 


k-1 


(2.6) 


where  p,  >  0  and  x,€[0,2tt).  Then  we  take  x,  ,  k  =  l,2,...,q  as  the  estimates 


of  x.  s. 
k 


Remark  2.1.  Sometimes,  u^+^  be  zero.  In  such  a  case,  there  may 
be  less  than  q  roots  for  B(z),  and  we  can  not  get  q  estimates  of  t^'s.  Howev 
in  the  large  sample  case,  we  can  prove  that  with  probability  cne,  u 


q+1 , q+1 


>  0 


for  large  N. 

Remark  2.2.  Using  the  Schmidt  orthogonalization  procedure,  we  can  seek 


for  0..  and  (u..  ) . 
N  lk 


6 


Remark  2.3.  Using  our  method,  we  do  not  bother  about  answering  the  two 
problems  mentioned  above. 

In  the  sequel,  we  will  establish  the  strong  consistency  of  x^'s.  Before 
doing  that,  we  introduce  the  following  lemma. 

Lemma  2.1.  Let  A  =  (a^)  anc^  ®  =  are  two  Hermitian  p*p  matrices  with 

spectrum  decompositions 


*  ** 

A  =  )  6.u.u.  ,  6,  >  6 >  ...  >6  , 
.L ,  i-i-i  1—  2  —  —  p 

i=l 


r  x 

B  =  )  A.v.v.,  A,  >  A_  >  ...  >  A  , 

.L  .  i-i-i  1  -  2  -  —  p 

1=1  r 

where  6's  and  A's  are  eigenvalues  of  A  and  B  respectively,  u's  and  v's  are 
orthogonal  unit  eigenvectors  associated  with  6's  and  A's  respectively.  Further, 
we  assume  that 


V  .+1  =  V  =  V  n0  =  0  "  nl  "  •••  *  n s  =  P.  h  =  1.2 . . 

h-1  h 


A,  >  A-  >  ...  >  A  , 
12  s 


and  that 


aik‘bikl  "  a*  i*  k  =  1 . P' 


Then  there  is  a  constant  M  independent  of  a,  such  that 

(i)  <  Mct*  1  *  1,2,.  ..,p 

(ii)  "h  u.u*  =  ft  v.v*  +  c(h)  with 
L  -1-1  L 

1=^.1+!  i-Vi+l 


c 


(h) 


(C 


(h) 

Hk 


),•  -|c 


(h) 

ilk 


£  Ma,'  2,',  k  =  1,2 


,  p ,  1)  1,2, ...  ,s. 


Proof.  By  Von-Neumann' s  inequality,  one  can  easily  obtain 

Y  (5. -A.)2  <  tr(A-B)2, 
i=l  1  1 


which  implies  (i)  with  M  =  p. 

For  simplicity,  we  denote  by  D  =  0(a)  the  fact  that  |d^l  <  Ma,  i=l,2,...,m, 
k  =  1,2,..., n  for  any  mxn  matrix  D  =  (d^).  To  prove  (ii),  without  loss  of 
generality,  we  can  assume 


s  - 


a  =  i  kI  vv  B  =  l  x  h  ^  w 

h-l  i6Ll  h-1  iei^-1"1 


where  =  {n^.j+l , . . . ,nh> .  When  s  =  1,  (ii)  is  trivial.  Now  we  assume  (ii) 
is  true  for  s  =  t-1,  and  proceed  to  prove  (ii)  for  s  =  t.  When  s  =  t, 
t-1  -  -  ...  t-1  -  - 


J/w.i  =  j/vv  i  vt +  oM- 

h=l  16L,  h=l  lfcL, 


(2.7) 


Multiply  from  right  hand  by  v^,  k€Lt  in  the  two  hand  sides  of  (2.7),  we  get 


which  implies  that 


t-1 


J<W  E  “i(“iV  =  0<01, 

h=l  16L, 


u.v, 

-l-k 


0(a),  i0Ls,  k6Ls> 


Thus,  we  have 


)V 

U1V2 


0(a),  VXU2  =  0(a), 


(2.8) 


Put  U2  =  V1G1  +  V2G2'  where  G1:  nt-lX^p_nt-l^’  G2:  ^p"nt-l^X^p-nt-l^‘  By  ^2‘8^ 


V2U2U2V2  =  V2(Ip-UlUt)V2  =  V2V2  +  0(a2) 


=  I  +  0(a) 

P"nt-1 


By  (2.8)  and  (2.9),  we  get 


0(a)  -  V*U2  -  Gx  +  V*V2G2  -  Gj, 


which  implies  that 


U2  *  V2G2  +  °(a)* 


By  (2.9)  and  (2.10), 


G2G2 


V2V2G2G2V2V2  *  V2U2U2V2  +  0(a) 


I  +  0(a) 

P“nt-1 


From  (2.10)  and  (2.11),  it  follows  that 


l  u  u  -  U-U  -  V  V  +  0(a) 
16L 


l  v  v  +  0(a), 
i6Lt  1 


(2.9 


(2.1 


(2.1 


(2.1 


and  that 
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t;2-  r*  *  _  - 

I  K  I  +  K.l  I 


h«l 


i6Lh 


i6Lt-l+Lt 


u.u. 

-1-1 


t-2 


I  A  I  v  v*  +  X  1  v  v  +  0(a) 

h-1  ni6Lh~1"1  c  1i6Lt_1+Lt'1~ 


By  the  induction  assumption, 


I  u.u  -  l  v  v  +  0(a),  h 
i6Lh  i6Lh 


1 , . . . ,t-2. 


(2.13) 


and 


i6L 


“i“i 


t-1 


i6Lt-l+Lt 


ViVi  + 


(2.14) 


Thus,  (ii)  is  true  for  s  *  t  by  (2.12)  -  (2.14).  Lemma  2.1  is  proved. 

We  have  the  following: 

THEOREM  2 . 1  Suppose  the  4*"  moments  of  s(t)  and  n(t)  are  finite.  Then 
the  estimates  t^'s  are  strongly  consistent. 


Proof.  Let  b  *  (b^ ,b2» . . • ,b^+^ ,0 . 0)  be  the  p*l  vector  whose  elements 


b^ ,b2, . . . »bq+i  are  the  coefficients  of  the  polynomial  bq+^  II  (z-e^Tk)  *  f(z) 


q+1 


k=l 


with  restrictions  I  lbkl  *  1  and  bq+1  >  0. 


k-1 


Let  -q+1  “  -»  Oq+2  *  (0,b1,b2,...,bq+1,0 . 0)  . Hp  -  (0 . 0,^,..., 


bq+^)^  be  all  p«l  vectors.  From  a^h^^  ■  ej<£"1)Tkf(ejTk)  -  0,  k  -  1,2,. ...q. 


t  ■  l,...,p-q,  it  follows  that  n  n  are  all  eigenvectors  of  ATA  associated 

-q+1  -p 


with  zero  eigenvalue.  Since  they  are  linearly  independent,  they  span  the  eigen- 


N 

subspace  of  A?A  associated  with  zero  eigenvalue.  Let  l ^  denote  this  subspace 


By  the 


and  let  ?2(u  +^)  =»  £  pP^n.  denote  the  projection  of  u 


k=q+l 

strong  law  of  large  numbers,  we  have 


,,  on  v_. 
q+1  2 


Z  -*■  I  **  A.'fk  +  o^I  ,  a.s.  as  N 
P 


by  Lemma  2.1,  it  is  easy  to  see  that 


Uq+i  =  p2^q+l^  +  a,s*  as  N 


(2. 


Since  Uq+j  ^  *  0  for  2.  =  q+2,...,p,  we  see  that  the  last  p-q-1  components 


of  ^2^q+l^  =  ^  tenc*  t0  zero  aljnost  surely.  From  this  and  the 


k*q+l 

expressions  of  n^*  k  =*  q+l,...,p,  noting  that  b^+^  are.  positive  constants, 
we  get 


(N) 

lim  p£  =0  a.s.  for  k  =  q+2,...,p 
N-h» 


and 


lim  P^?  *  1  a.s., 

N~  q+1 


which  implies  that 


u  . ,  -*■  h  . .  =  b,  a.s.  as  N  +  ». 
-q+1  -q+1 


(2. 


By  the  definition  of  b,  we  know  that  e^Tk,  k  *  l,2,...,q,  are  the  roots 
of  the  polynomial  equation 

1±1 


Y- 


k«l 


bfcz  x  *  0. 


(2.] 


Hence,  after  suitable  rearrangement 
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and  consequently, 


-  JTk  JTk 

p^e  ■+  e  ,  a.s.,  k  *  l,2,...,q. 


Pk  •*  1  a.s. 


■*  k  =  l,2,...,q,  a.s.  as  N  ■* 


(2.17) 


which  proves  the  theorem. 


Remark  2.4.  If  q  is  known,  to  ensure  the  strong  consistency  of  x^'s, 
we  only  need  to  assume  the  second  moments  of  s(t)*s  and  n(t)'s.  But  in  ITC 
procedure,  to  guarantee  the  strong  consistency  of  the  estimate  of  the  signal 
number,  we  assumed  the  A1"*1  moments  of  s(t)  and  n(t)  exist.  (Refer  to  Zhao, 
Krishnaiah  and  Bai  [1986a]).  Therefore  in  this  theorem  we  still  assume  the 
4fck  moments  exist,  so  that  the  conclusion  of  Theorem  2.1  is  still  true  by 
using  the  ITC  estimate  of  signal  number  q  instead  of  q  when  q  is  unknown. 
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3.  ESTIMATE  OF  DOA  IN  THE  PRESENCE  OF  COLORED  NOISE 

In  the  section  2,  we  obtain  an  estimate  of  DOA's  when  the  additive  sensor 
noise  is  spatially  white.  When  the  sensor  noise  is  colored,  the  case  is  more 
complicated.  For  this  case,  Paulraj  and  Kailath  (1986)  •  proposed  a 
solution  to  the  DOA  estimation  problem.  Their  technique  is  applicable  to 
situations  where  it  is  possible  to  obtain  two  estimates  of  the  array  covariance 
in  which  the  unknown  noise  field  remains  invariant  while  the  signal  field  under¬ 
goes  some  change.  This  method  is  based  on  computing  the  difference  of  the  two 
measured  covariances,  thus  subtracting  out  the  unknown  noise  covariance  and 

leaving  only  the  difference  matrix  of  the  two  signal  covariance. 

Assume  that  there  are  two  estimates  of  the  array  covariance  with  the 

array  being  displaced  between  the  measurements.  This  displacement  could  be 

of  several  types.  Examples  of  spatial  displacements  are  rotations,  translations, 

or  a  combination  of  the  two.  Displacements  can  also  be  of  a  temporal  nature 

with  the  noise  statistics  being  long-term  stationary,  while  those  of  the  signals 

are  only  short-term  stationary.  For  the  details,  refer  to  Paulraj  and  Kailath 

(1986).  Here  we  assume  that  the  noise  covariance  matrix  is  invariant  across 

the  two  measurements  while  the  signal  covariance  matrix  and  the  DOA.' s  change  in 

some  manner  between  the  measurements.  Thus  we  have  the  following  model: 

x(g,)(t)  *  A(£)s(8,)(t)  +  n(i)(t),  t  =  1,2,. ..,N,  i  =  1,2  (3.1) 


where  x^^(t):  p*l,  the  observations  received  by  p  sensors  for  the  4*"*1  measurements. 


(4) 

s  (t):q^*l,  the  signal  vector  emitted  by  q^  sources. 


1,2,  n^\t)  is  the 


colored  noise  for  the  4fc^  measurements, 


( 1  ,e”jT^,  ^ .  .  ,e'^p"1)Tk  ))T,  k  -  1,2 . q:  and  4 


(a{l\...,a(^)  and  a££)  ' 
-1  -qj  -k 


1,2. 


r. 


It  is  usual  to  assume  that 

(i)  For  each  £,  £  =  1,2,  {sv  '(t)}  iid.,  {nv  '(t)}  iid.,  and  independent 
of  (sU)(t)}, 

(ii)  Es(£)(t)  -0,  En(£)(t>  =  0,  / 

Es(£)(t)s(i)*(t)  =«  >0,  /  (3.2 

En(£)(t)n(£)*(t)  *  I,0  >  0,  £  =  1,2,  t  =  1,2 . N,  l 

where  and  z ■  are  all  unknown. 

(£ ) 

The  covariance  matrix  of  x  (t)  is  given  by 


2-U)  «  A.U)f(*)A(0*  +  z  m  1>2. 

For  the  translational  invariance  model,  we  know  that  *  A^2\  and 


Z(D  _  j.(2)  ,  Ad >(^( 1  )_^(2 ) )A( 1  >* , 


where  we  assume  that  ^g  Qf  ranjl[  ancj  <  p_  Also,  we  assume 

that  rf^'s  are  distinguished.  This  means  that  A^  is  of  ful-1  -rank  (i.e., 


For  other  invariance  models,  we  have 


?(D  0 


Z(1)  -  I(2)  «  (A<l>.A<»^  _y(2))(A(1),A(2))*, 


were  we  assume  that  x^  's  and  x^  's  are  all  distinguished  and  2q^  <  p.  In 
this  case,  (A^^,A^^)  is  full  rank  (i.e.,  ■  2q^). 

For  the  model  (3.3),  we  write  A  *  A^\  T  ■  -  4^^,  and  q  *  q^. 

For  the  model  (3.4),  we  write  A  *  (A^*\a^),  T  ■  ®C2^  and  q  ■  2q^. 


Et  -  £  l  xU)(i)xU)*(i),  £  -  1,2. 

i-1 


1A 


Let  A.,  >  A._  >  . . .  >  A.  >  A  A  *0  and  6.  >  6_  >  ...  >  6  denote 

1  —  2  -  -  q  q+1  p  1  ~  2  —  -  p 


2  »  ^  2 

the  eigenvalues  of  (Z^Z^  and  (Z1“Z2)  respectively.  Take  CN  satisfying 


CN  n  CN  . 

TT*0  and  loglogii  as  M-"- 


(3.5) 


Write 


I(k,C  )  -  N  l  6  +  k  C  , 

”  i-k+1 


(3.6) 


and  define  q  as  follows 


l(q,CN)  -  min{I(0,CN) . I(p-1,CN)}. 


(3.7) 


We  have  the  following 


THEOREM  3.1.  Suppose  that  (3.2)  holds,  A  is  of  rank  q  and  the  A  moments 


of  n^  ^(t)  exist  for  l  m  1,2.  Then  q  is  a  strongly  consistent  estimate  of  q. 


Proof.  By  the  law  of  the  iterated  logarithm,  for  l  ■  1,2, 


5^  +  0('//^loglogN)  a.s.  as  N  ■»  », 


(3.8) 


Us ing  Lemma  2.1,  we  have 


lim  6,  *  A,  a.s.,  i  *  1,2,. ...p. 
N-*» 


(3.9) 


Since  the  matrix  Z^-  Is  of  rank  q,  there  exists  a  p*(p-q)  matrix  QqQq  *  I  _£ 


and  Qq(Z^-Z2)  *  0*  It  is  well  known  that 


l  6  -  min  trQT(i . -£7)2Q. 

i-q+1  JO-1*  a 

p-q 


(3.10) 


By  (3.8), 


!y 

v 


& 


fe: 


K 


H 


Qq(z1-I2)  *  Qq  (  zi-Z2"^5‘1~5:2^  “  OC/iloglogN).  a.s.  (3.11) 

By  (3.10)  and  (3.11), 

0  1  PI  1  tr  Qj(  Z  -L)2Qq  *  O(iloglogN),  a.s.  (3.12) 

i*q+l 

Using  (3.9)  and  (3.12),  noticing  that  A^  >  0,  we  can  easily  prove  that,  with 
probability  one  for  N  large, 

l(q,CN)  <  I(k,CN),  k  *  q,  k  <  p-1, 

which  implies  that 

q  -  q. 

Theorem  3.1  is  proved. 

-  *  -  2 

In  the  sequel,  we  assume  that  q  is  known.  Write  T  =  TA  AT,  then  (l^-Z^) 
can  be  rewritten  as 

(vy2  *  ata*. 

Note  that  A  is  of  the  form  A  ■  (a,,..., a  )  with 

-1  -q 

-jTk  -j(p-lhk  T 

^k  *  ( 1 »6  » • • »  >6  )  t  ^  1 » • • • 

where  x^'s  are  distinguished.  So  the  problem  of  estimating  the  DOA's  reduces 

«  «  2 

the  case  of  section  2.  Let  u^,...,Up  denote  the  eigenvectors  of  (Z^-S^) 
associated  with  6^,...,<5p.  Based  on  W^  *  (u^+^ , . . . ,u^) ,  we  can  use  the  method 
proposed  by  us  in  the  section  2,  and  take  x^,  k  -  1,2 
x^'s.  In  the  same  way,  we  have 


, ...,q,  as  the  estimates  of 


of  Theorem  3.1,  t^'s  are  strongly 


unclassified 


mm 


— i 
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